Dramatic reductions in inelastic cross sections 
for ultracold collisions near Feshbach resonances 
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We show that low-energy inelastic cross sections can decrease as well as increase in the vicinity of 
a zero-energy Feshbach resonance. When an external field is used to tune across such a resonance, 
the real and imaginary parts of the scattering length show asymmetric oscillations, with both peaks 
and troughs. In favorable circumstances the inelastic collision rate can be reduced almost to zero. 
This may be important for efforts to achieve evaporative and sympathetic cooling for molecules. 



Cold and ultracold molecules have fascinating proper- 
ties that will find applications in many areas of physics, 
ranging from precision measurement to quantum com- 
puting [l| . Ultracold molecules offer new possibilities for 
quantum simulations and quantum control, while quan- 
tum gases of ultracold polar molecules are expected to 
exhibit a wide range of new quantum phases. 

Cold and ultracold molecules must always be confined 
in traps, and trap losses are crucial. In particular, coUi- 
sional stability is very important. Magnetic and electro- 
static traps can trap molecules only when they are in low- 
field-seeking states, and such states are never the lowest 
in the applied field. Any inelastic collision that transfers 
internal energy into relative translational motion causes 
either heating or trap loss. It is thus very important to 
understand inelastic collisions and to find ways to min- 
imise them. The purpose of the present paper is to show 
that inelastic collision rates can sometimes be dramati- 
cally reduced by tuning close to a Feshbach resonance [3| 
with an applied electric or magnetic field. 

It has been possible for about 5 years to create 
molecules in highly vibrationally excited states in ultra- 
cold atomic gases [Sj, both by photoassociation and 
by magnetoassociation [a]. A major goal was achieved 
in 2008 when Ni et al. [y] succeeded in transfering KRb 
molecules formed by magnetoassociation at T 350 nK 
into their ground rovibrational state by stimulated Ra- 
man adiabatic passage (STIRAP). Experiments on coUi- 
sional trap loss are already under way 01 . Danzl et al. 
[1,0] and Lang et al. [10] have carried out analogous ex- 
periments on Cs2 and triplet Rb2, respectively. There 
have also been considerable successes in direct photoas- 
sociation to produce low- lying states [ll|j d) [U [3] • 

Methods based on photoassociation and magnetoasso- 
ciation are limited to molecules formed from atoms that 
can be laser-cooled, such as the alkali metals. However, 
a wider range of molecules can be cooled directly from 
high temperature to the millikelvin regime, using meth- 
ods such as buffer-gas cooling and Stark deceleration 
. Polar molecules such as ND3 and OH can be held in 
electrostatic or alternating current traps [3 ' while 
paramagnetic molecules such as CaH, O2, NH and OH 



can be held in magnetic traps. However, at present the 
lowest temperatures that can be achieved in static traps 
are around 10 mK, and there are a variety of propos- 
als for ways to cool them further, including evaporative 
cooling, sympathetic cooling and cavity-assisted cooling 

Magnetic fields have important effects on the interac- 
tions and collisions of paramagnetic molecules [13, [H, 

. In previous work [2^ , we explored the use of Fesh- 
bach resonances to control molecular collisions with ap- 
plied fields. For the prototype system He -I- NH ('^S"), 
we located magnetic fields at which bound states cross 
open-channel thresholds and then characterized the re- 
sulting low-energy Feshbach resonances as a function of 
magnetic field. For a resonance at which a bound state 
crosses the lowest open-channel threshold, the real scat- 
tering length a behaves in the same way as in the atomic 
case [26] and exhibits a pole as a function of magnetic 
field B. However, for resonances in which a state crosses 
a higher threshold, we observed quite different behavior. 
In this case, inelastic scattering can occur and the scat- 
tering length is complex. The complex scattering length 
a{B) = a{B) — i(3{B) was found to follow the formula 



a[B) = flbg - 



2(B-B,cs)/r'^"' + i 



(1) 



where Ohg is a slowly varying background term and 
Bros and F'^"' are the position and width of the reso- 
nance. The resonant scattering length 0^3 characterizes 
the strength of the resonance. The elastic and total in- 
elastic cross sections are given approximately by 



aei(B) «4^|a(B)|2 and al°l,{B) 



47r/3(B) 
fco 



(2) 



The inelastic scattering in He + NH is very weak ex- 
cept near resonance. Under these circumstances, aci{B) 
shows a symmetric oscillation at resonance and al°l^{B) 
shows a simple peak [2^ . This corresponds to a real value 
of flics- However, a more complete derivation [27] subse- 
quently showed that, when the background scattering is 
significantly inelastic, it is possible for flros to be complex. 



2 



When this occurs, inelastic cross sections show troughs 
as weh as peaks near resonance. This is potentially of 
great importance, since inelastic collisions generally pro- 
vide trap loss mechanisms and strong inelastic processes 
can prevent evaporative or sympathetic cooling. 

In this letter we consider a more strongly coupled 
system, with significant inelastic scattering, in order to 
demonstrate the dramatic reductions in inelasticity that 
can occur near Feshbach resonances. The system we have 
chosen is ^He + "'^^02 for which a reliable poten- 

tial energy surface has been calculated by Groenenboom 
and Struniewicz j2^. The ^^©2 molecule has a ground 
state with rotational quantum number n = 1, because 
the ^^O nucleus is a boson with nuclear spin / = 0. 
The bound-state Schrodinger equation for ''He -I- ^^©2 is 
solved by propagating coupled differential equations us- 
ing the BOUND package f29l . Isoj , as modified to handle 
magnetic fields [25|. The calculations are carried out in 
a completely decoupled basis set [2^, \nmn)\sms}\LmL) , 
where s = I is the electron spin of O2 and L is the end- 
over-end rotational angular momentum of He and the 
molecule. All the m quantum numbers represent space- 
fixed projections on the axis defined by the magnetic 
field. The only good quantum numbers are the par- 
ity (—1)"+^+! and the total projection quantum num- 
ber Mtot = rrin + rris + ttil- At energies above the low- 
est threshold, BOUND locates both physical quasibound 
states and artificial states that result from box-quantizing 
the continuum for open channels. However, it is straight- 
forward to identify the physical states by inspecting the 
dependence of the eigenvalues on the outer limit of the 
propagation. 

Fig. [1] shows the bound and quasibound states of the 
^He-^^02 complex near the n = 1 thresholds as a func- 
tion of magnetic field B, with artificial levels removed, 
together with the thresholds for dissociation to form He 
+ O2. The thresholds are characterized at zero field by 
O2 quantum numbers tt., s,j with s = I and j = 0, 1,2 
for n = 1, and each one splits into 2j + 1 components 
labeled by mj at non-zero magnetic field. He-02 is a 
weakly anisotropic system, so n and s remain nearly 
good quantum numbers and the levels of the complex 
are characterized by additional quantum numbers L and 
J, where the total angular momentum J is the resul- 
tant of j and L. At zero field each {n,j,L) level splits 
into min(2j -|- 1, 2L -t- 1) sublevels with different values of 
J. When a magnetic field is applied, each sublevel splits 
into 2J+1 components with different values of A/tot- The 
J quantum number remains a useful label for magnetic 
fields up to about 1000 G, but above that the levels of 
different J are strongly mixed. By about 5000 G the lev- 
els have separated into groups that may be labeled with 
an approximate quantum number fhj that takes values 
from —j to +j. 

Crossings between quasibound states and thresholds 
will produce zero-energy Feshbach resonances in s-wave 
scattering if an L = scattering channel is permitted by 
the constraints on parity and Mtot- This occurs only for 
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FIG. 1: (Color online). The pattern of levels from bound-state 
calculations on ^He-^''02 near the n = 1 thresholds, with ar- 
tificial levels removed, as a function of magnetic field B. The 
calculations are for even parity, Mtot = -6 to +6. The "O2 
threshold energies are shown as solid black lines. The cir- 
cles show crossings between bound states and thresholds with 
nij = Mtot that produce zero-energy Feshbach resonances in 
s-wave scattering. 



thresholds corresponding to mj = Mtot as shown by the 
circles in Fig. [TJ In the present work we are particularly 
interested in resonances that occur at eiczied thresholds, 
where inelastic scattering may occur. 

Once the crossing points have been located in Fig. [U 
we carry out scattering calculations, holding the kinetic 
energy fixed at a small value while sweeping the mag- 
netic field across the resonance. This is done using the 
MOLSCAT packa ge [ 3ll. as modified to handle collisions 
in magnetic fields [25j . 

Typical resonance profiles for *He -I- O2 are shown in 
Fig. [2I for the resonance labeled 1 in Fig. [T] at collision 
energies E ^ 1 /iK, 100 fxK and 10 mK. At all three en- 
ergies the total inelastic cross section (summed over out- 
going partial waves L') drops by almost a factor of 1000 
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FIG. 2: (Color online). Elastic (red) and total inelastic 
(green) cross sections for resonance 1 in He + O2 as a function 
of magnetic field at collision energies of 10"^ K (solid lines), 
10""* K (dashed lines) and 10"^ K (dotted lines). 
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FIG. 3: (Color online). Elastic (red) and total inelastic 
(green) cross sections for resonance 2 in He + O2 as a function 
of magnetic field at collision energy 10~^ K. 



from its background value at a field just below Bros- At 
higher energies the resonance is shifted slightly, but the 
resonant suppression is just as strong. The resonant con- 
tribution to the inelastic cross section follows the 
Wigner threshold law [3^ only at very low collision en- 
ergies (below 100 /iK). At higher energies it actually de- 
creases faster than predicted by the threshold law and the 
elastic cross section also decreases. The p-wave contri- 
bution to the inelastic cross section is non-resonant, and 
increases approximately as E'^^^^ with energy while the 
s-wave contribution decreases. For the resonance in Fig. 
[2] its value is about 1 at 10 mK, so that p-wave scat- 
tering will dampen the suppression of inelastic scattering 
at temperatures above this. 

The resonances for He -I- O2 are quite wide, with 
IFbI = 250 to 500 G. The resonance in Fig. [U pro- 



vides substantial suppression of inelastic cross sections 
across a range of at least 100 G. This would be suffi- 
cient to provide a working energy range of about 10 mK 
for magnetically trapped states of a molecule in a 
electronic state. The broad resonances observed here 
contrast with those previously characterized for He-NH, 
which had \Tb\ < 10~^ G [2l|. The difference arises 
in this case because the n = 1 closed channels involved 
for O2 are directly coupled to the inelastic channel(s) 
by the potential anisotropy, whereas the n = closed 
channels involved for NH were only indirectly coupled to 
the open channels. However, broad resonances are likely 
to be common in molecule-molecule systems, and indeed 
in atom-molecule systems involving atoms heavier than 
He, because the couplings due to potential anisotropy are 
much stronger in heavier systems [33l [33 | . 

The asymmetric lineshapes observed here are analo- 
gous to Fano lineshapes [35] in bound-free absorption 
spectra. Fano considered the interference between the 
bound and continuum contributions to a transition ma- 
trix element near resonance. He showed that the bound- 
state contribution rises from zero to a peak at reso- 
nance while the continuum contribution drops from its 
background value to zero and changes sign at resonance. 
When there is only a single continuum channel, there is 
always a point near resonance where the bound and con- 
tinuum contributions cancel completely. However, when 
there are N outgoing channels, there is one particular 
linear combination of them that is coupled to the bound 
state and A'^ — 1 orthogonal linear combinations that are 
not [3^. The resonance suppresses inelastic scattering 
into the former but not into the latter, so the cross sec- 
tion does not drop to zero. 

For low-energy resonant scattering in the presence of 
inelastic channels, the partial width for the incoming 
(elastic) channel is proportional to the incoming wavevec- 
tor fco, while the partial widths for the inelastic channels 
are essentially independent of fco 27[. At low energies 



we may therefore consider the bound state to be coupled 
only to the outgoing (inelastic) channels and apply Fano 
theory directly to the inelastic cross sections. 

Even for ^^©2 molecules at the j — 2, nij = — 2 thresh- 
old, which can relax only to form j = 0, rrij = 0, there are 
outgoing channels with several values of L' . For s-wave 
scattering [L = 0), L' must be at least Ato^ and must 
be even to conserve parity. The kinetic energy release of 
1.9 K at i? = 9750 G is above the centrifugal barrier for 
L = 2 (0.4 K) but below that for L = 4 (2.4 K). Because 
of this, the L — 2 channel dominates the inelastic scat- 
tering away from resonance and is the outgoing channel 
most strongly coupled to the bound state. The inelastic 
cross section therefore shows a deep minimum, though 
there is still a little background inelastic scattering that 
is not suppressed by the resonance. 

The situation is somewhat different for resonances at 
the j = 2,mj = — 1 threshold, such as that shown in 
Fig. [3] (resonance 2 in Fig. [T]). In this case the res- 
onance suppresses the total inelastic cross section by 
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less than a factor of 10 from its background value. At 
11660 G the kinetic energy release is 0.55 K for relax- 
ation to the j — 2, rrij = —2 (upper) threshold and 2.3 
K for relaxation to the j — 0, rrij — (lower) thresh- 
old. The L — 2 outgoing channels at both inelastic 
thresholds contribute significantly to the inelastic scat- 
tering far from resonance. However, the resonant bound 
state (with rrij — +2) is coupled much more strongly 
to j = 0, mj = channels than to j = 2, ruj — —2 
channels. The resonance therefore suppresses inelastic 
scattering into the lower channel but there is significant 
background inelastic scattering into the upper channel 
that is unaffected by the resonance. 

The kinetic energy release needed to surmount cen- 
trifugal barriers depends on the reduced mass and will 
be smaller in heavier systems than for He -I- O2, as will 
the temperature range in which s-wave scattering dom- 
inates. However, much lower temperatures are already 
being achieved in current experiments @, Q . 
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this will be systems in which all but one of the outgoing 
channels are suppressed by centrifugal barriers. 
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